Universal correlations and coherence in quasi-two-dimensional trapped Bose gases 
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We study the quasi-two-dimensional Bose gas in harmonic traps at temperatures above the 
Kosterlitz-Thouless transition, where the gas is in the normal phase. We show that mean-field 
theory takes into account the dominant interaction effects for experimentally relevant trap geome- 
tries. Comparing with Quantum Monte Carlo calculations, we quantify the onset of the fluctuation 
regime, where correlations beyond mean-field become important. Although the density profile de- 
pends on the microscopic parameters of the system, we show that the correlation density (the 
difference between the exact and the mean-field density) is accurately described by a universal ex- 
pression, obtained from classical-field calculations of the homogeneous strictly two-dimensional gas. 
Deviations from universality, due to the finite value of the interaction or to the trap geometry, are 
shown to be small for current experiments. We further study coherence and pair correlations on a 
microscopic scale. Finite-size effects in the off-diagonal density matrix allows us to characterize the 
cross-over from Kosterlitz-Thouless to Bose-Einstein behavior for small particle numbers. Bose- 
Einstein condensation occurs below a characteristic number of particles which rapidly diverges with 
vanishing interactions. 

PACS numbers: 03.75.Hh,05.30. Jp 



I. INTRODUCTION 

In recent years, several experiments [1, 2] studied two- 
dimensional ultra-cold atomic gases from the normal 
phase down in temperature to the Kosterlitz-Thouless 
transition [3] and into the low-temperature superfluid 
phase. The interference of two simultaneously prepared 
two-dimensional gases evidenced the presence of vortices 
[1]. Related experiments investigated interaction and 
correlation effects [2, 4] in the density profile and in co- 
herence patterns. For a quantitative description of the 
Kosterlitz-Thouless transition and of the interaction ef- 
fects, it proved necessary to account for the quasi-two- 
dimensional nature of the gas, that is to include thermal 
excitations in the strongly confined z-axis in addition to 
the weak trapping potential in the xy-plane [5, 6]. 

In weakly interacting two-dimensional Bose gases, 
the Kosterlitz-Thouless phase transition occurs at rel- 
atively high phase-space density (number of atoms per 
phase-space cell A T = 2irh 2 /mT). This density is 
n c A T ~ log(£„/<?) [7-9], where g characterizes the two- 
dimensional interaction strength, T is the temperature, 
m the mass of the atoms, and n the density. The coef- 
ficient £ n = 380 ± 3 was determined numerically using 
classical-field simulations [9]. For the ENS experiment 
of Hadzibabic et al. [1], the critical phase-space den- 
sity is n c A T r~j 8 in the center of the trap, whereas in 
the NIST experiment of Clade et al. [2], n c A T is close 
to 10. For phase-space densities between one and the 
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critical number, the gas is quantum degenerate yet nor- 
mal. The atoms within one phase-space cell are indistin- 
guishable. They lose their particle properties and acquire 
the characteristics of a field. The mean-field description 
of particles interacting with a local atomic density n(r) 
may further be modified through correlations and fluctu- 
ations. Quantum correlations can be several times larger 
than the scale At- This gives rise to "quasi-condensate" 
behavior inside the normal phase. 

In this paper, we study the quantum-degenerate regime 
at high phase-space density in the normal phase. We 
first discuss the peculiar quasi-two-dimensional thermo- 
dynamic limit where, as the number of particles in the 
gas is increased, the interactions and the lattice geome- 
try are scaled such that a finite fraction of all particles 
are in the excited states of the system. In this thermo- 
dynamic limit, the Kosterlitz-Thouless transition takes 
place at a temperature comparable to the Bose-Einstein 
transition temperature in the non-interacting case, and 
the local-density approximation becomes exact. We first 
clarify the relation between different recent versions of 
quasi-two-dimensional mean-field theory [6, 10, 11] in the 
local-density approximation (LDA), and also determine 
the finite-size corrections to the LDA. We compare mean- 
field theory to a numerically exact solution obtained by 
path-integral Quantum Monte Carlo (QMC) calculations 
with up to N > 10 5 interacting particles in a harmonic 
trap with parameters chosen to fit the experiments. We 
concentrate on the correlation density, the difference be- 
tween the exact density and the mean-field density at 
equal chemical potential, and show that it is essentially 
a universal function, independent of microscopic details. 
Within classical-field theory, the correlation density is 
obtained from a reparamctrization of known results for 
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the strictly two-dimensional homogeneous system [12]. 
The classical-field results hold for small interaction pa- 
rameters g — ► 0, but our full QMC solution accounts 
for corrections. We compute the correlation density by 
QMC and show that it is largely independent of the trap 
geometry, the temperature, and the interaction strength. 

We also study off-diagonal coherence properties, and 
the density-density correlation function of the quasi-two- 
dimensional gas. It is well known that even at high 
temperature, bosonic bunching effects enhance the pair- 
correlation function on length scales below At which for 
the ideal Bose gas approaches the characteristic value 
2n 2 at vanishing separation. In our case, interference in 
the z-direction reduces the in-plane density fluctuations 
even for an ideal gas and within mean-field theory, and 
the reduction of the pair correlations from 2n 2 no longer 
proves the presence of beyond-mean-field effects. 

We finally discuss finite-size effects in the quasi-two- 
dimensional Bose gas. For the density profile, they are 
not very large, but we point out their great role for 
off-diagonal correlations. The latter are responsible for 
a cross-over between the physics of Bose-Einstein con- 
densation at small particle number and the Kosterlitz- 
Thoulcss physics for larger systems; both regimes are of 
relevance for current experiments. This cross-over takes 
place at a particle number TV 
rapidly as the interaction in the gas diminishes 
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II. SYSTEM PARAMETERS AND MEAN-FIELD 
DESCRIPTION 

A. Quasi-two-dimensional thermodynamics 

We consider TV bosons in a three-dimensional pancake- 
shaped harmonic potential with parameters lo x = uj v = lu 
and io z 3> lu at inverse temperature (5 = 1/T. The z 
variable is separate from x and y, and wc denote the 
three-dimensional vectors as r — (r, z), and write two- 
dimensional vectors as r = (x,y), and r = |r|. 

The quasi-two-dimensional regime of the Bose gas [6] is 
defined through a particular thermodynamic limit TV — ► 

00, where the temperature is a fixed fraction t = T/Tj^ c 
of the Bose-Einstein transition temperature of the ideal 
two-dimensional Bose gas, T|g C = V '6TV 'Hlu / 'it . Although 
the two-dimensional trapped Bose gas undergoes a Bose- 
Einstein transition only for zero interactions, still 
sets the scale for the Kosterlitz-Thouless transition in 
the interacting gas[6, 8]. In the quasi-two-dimensional 
regime, a finite fraction of atoms remains in excited states 
in z. The excitation energy is scaled as Hlu z oc T^ c , 
which implies that lu z increases as TV 1 / 2 in the thermo- 
dynamic limit. 

Interatomic collisions are intrinsically three- 
dimensional. Here, we consider the experimentally 
relevant case where the range of the scattering potential 
ro is much smaller than the typical lateral extension 

1, = {musz/h) -1 / 2 , and also where rg is much smaller 



than the inter-particle distances. The interactions 
are then described by the three-dimensional s-wave 
scattering length a s , and one may characterize the 
quasi-two-dimensional gas through a bare effective 
two-dimensional interaction strength, g, 



4nh 2 a s 

g = , 



(1) 

(2) 



where ipo( z ) is the unperturbed ground state of the con- 
fining potential and g is the usual three-dimensional 
coupling constant. For a harmonic confinement, g = 
\/87ra s /Z z , and a s /l z must be kept constant in the ther- 
modynamic limit to obtain a fixed two-dimensional in- 
teraction strength. 

Quasi-two-dimensional scattering amplitudes depend 
logarithmically on energy, e, in terms of a universal func- 
tion of a s /l z and of e/huj z ~ T/u> z , which are both kept 
constant in the quasi-two-dimensional thermodynamic 
limit. The logarithmic energy dependence yields small 
corrections of order (a s /l z ) 2 [13-15] to the bare interac- 
tion g. They can be neglected in the following. 

The scaling behavior in the quasi-two-dimensional 
limit corresponds to the following reduced variables: 



r/l T z = z/l z 



t = r/T, 



2d 
BEC 



nw z /T 2d 



BEC 



(3) 



g = mg/{V2Trl z H 2 ) 



where It = (T/muj 2 ) 1 / 2 is the thermal extension in the 
plane. The quasi-two-dimensional limit consists in tak- 
ing TV — ► oo, with t, uj z and g all constant. In this limit, 
Zt/At = t\/3TV/7r 3 ^> 1, so that macroscopic and mi- 
croscopic length scales separate, and the scaling of the 
three-dimensional density, 713^, is at constant 

n 3 d = md^Tlz- (4) 

In reduced variables, the normalization condition TV = 
J dr n(r) = (Zt/At) 2 / dr n(f) is expressed as 



df rn(f) 



m 2 ' 



(5) 



The local-density approximation becomes exact in the 
quasi- two-dimensional limit. 

In the ENS experiment, 87 Rb atoms are trapped at 
temperatures r ss 50 - 100 nK. The in-plane trap- 
ping frequencies are lu / (2n) ~ 50Hz whereas the con- 
finement is of order lo z /{2it) w 3kHz. With TV - 2 • 10 4 
atoms trapped inside one plane, typical parameters are 
t bec ~ 300nK (using h/k B ^ 7.64 • 10" 3 nKs), so that 
Cj z f=a 0.44 — 0.55. The scattering length a s = 5.2nm 
leads to an effective coupling constant g = 0.13, using 
h/m ~ 6.3 • 10~ 8 m 2 s~ 1 A~ 1 , where A is the atomic mass 
number. In the NIST experiment, sodium atoms at T ss 



2d 

BEC 



100nK arc confined by harmonic trapping poten- 
tials with w/(2tt) w 20Hz, and uj z /(2tt) ks 1kHz. This is 
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described by reduced parameters g = 0.02 and Q z = 0.50. 
The critical densities are n c sa log(380/<?) « 8.2 for the 
ENS parameters, somewhat lower than the NIST value 
h c w 9.9. Using the quasi-two-dimensional mean-field es- 
timates of Ref. [6] , the Kosterlitz-Thouless temperatures 
are located at <kt = Tkt/T^ c « 0.69 and £kt ~ 0.74, 
respectively. 



The quasi-two-dimensional limit describes a kincmat- 
ically two-dimensional gas, whose extension in the z di- 
rection is of the order of the thermal wavelength At- As 
u) z = \\ j {2-Ktl 2 z ) is decreased, a system at finite N turns 
three-dimensional. This is already the case for the ideal 
quasi-two-dimensional gas (with g = 0) where the Bosc- 
Einstein transition temperature crosses over from two- 
dimensional to three-dimensional behavior as a function 
of Q z , with asymptotic behavior given by 



l BEC 



C(2) 
W) 



1/3 



-1/3 



1C(2)~ 



6 03) 
2^2)3/2 CX P (-W z ) 



for G) z < 1 
for (D z 3> 1 



(6) 



(see [6]). In Eq. (6), the first term for ui z <C 1 de- 
scribes three-dimensional Bose-Einstein condensation in 
an anisotropic trapping potential. 

For the interacting Bose gas, the nature of the 
Kosterlitz-Thouless transition in two dimensions dif- 
fers from the Bose-Einstein transition of the three- 
dimensional gas. For large u z , universal features of 
the Kosterlitz-Thouless transition are preserved, but the 
density profiles and the value of the Kosterlitz-Thouless 
transition temperature depend on Cj z and g [5, 6]. For 
small confinement strength uj z , a dimensional cross-over 
between the two-dimensional Kosterlitz-Thouless transi- 
tion and the three-dimensional Bose-Einstein condensa- 
tion takes place at particle numbers such that the level 
spacing in the confined direction is comparable to the 
(two-dimensional) correlation energies, tuj z < gn/ir . 

The quasi-two-dimensional limit differs from the "ex- 
perimentalist's" thermodynamic limit where the atom 
number is increased in a fixed trap geometry, and at 
constant temperature. In this situation, the ratio be- 
tween the microscopic and the macroscopic length scales, 
It/Xt = T/(tkJy/2n), remains constant and finite. The 
number of particles in any region of nearly constant den- 
sity remains also finite so that, in contrast to the quasi- 
two-dimensional thermodynamic limit, corrections to the 
LDA persist. 



B. iV-body and mean-field Hamiltonians 

The gas specified in Section II A is described by the 
Hamiltonian 



H =H + V, 

N 



i=l 
N 

V = v(\ri-r. 
%<j=i 



^ 2 Vf 1 r 2 2 



2^21 



(7) 

(8) 

(9) 



where v is the three-dimensional interaction potential. 
We compute the iV-body density matrix at finite temper- 
ature using three-dimensional path-integral QMC meth- 
ods. We thus obtain all the thermodynamic observ- 
ables [5, 16, 17] for up to N = 10 6 . QMC calculations 
have clearly demonstrated the presence of a Kosterlitz- 
Thouless transition [6] for parameters corresponding to 
the ENS experiment. 

In the mean-field approximation, one replaces the N- 
body interaction between atoms in Eq. (9) by an effec- 
tive single-particle potential. The mean-field Hamilto- 
nian writes 



H 



mf 



Hn + V, 



mf 5 



(10) 



where 



V, 



N 



mf 



2gn 3d (f i) -g df [n 3d {P)} 2 (11) 



is the mean-field potential energy. From the corre- 
sponding partition function in the canonical or grand- 
canonical ensemble, all thermodynamic quantities can 
be calculated. The three-dimensional density n 3 d(P) in- 
side the mean-field potential must be determined self- 
consistcntly. In all situations treated in the present pa- 
per, self-consistency is reached through straightforward 
iteration. 

Mean-field theory leads to an effective Schrodinger 
equation for the single-particle wavefunction, ^ (r ), of 
energy e,-, 



2yy2 



1 



2 2 

mu) r 



tfd 2 1 2 2 „ 

2m ' 2 2^ + 2 m ^ Z 

*,•(?) = e^f) (12) 

together with the total density 

m{p)*j(f) 



n 3 d(r) = J2 



(13) 



The exact solution of the mean-field eigenfunctions and 
eigenvalues for finite systems, is rather involved, but con- 
siderably simplifies in the local-density approximation. 

At finite N, in the canonical ensemble, we solve the 
mean-field equations through a Quantum Monte Carlo 
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simulation with N particles which avoids an explicit cal- 
culation of all eigenfunctions. In contrast to the usual 
interaction energy within QMC, consisting, in general, 
of a pair interaction potential, the mean-field interaction 
energy is simply given in terms of an anisotropic, single- 
particle potential proportional to the three-dimensional 
density profile, n3d(f), as in Eqs (10) and (11). This 
interaction potential must be obtained self-consistcntly 
as usual in mean-field. Once self-consistency in the den- 
sity is reached, one can compute correlation functions 
and off-diagonal elements of the reduced one-body den- 
sity matrix. 



C. Mean-field: Local-density approximation 

Mean-field theory simplifies in the quasi-two- 
dimensional thermodynamic limit, as the local-density 
approximation then becomes exact. This is because 
the natural length scale of the system, At, separates 
from the macroscopic scale It of variation of the density 
(Xt/It — > 0). The particle numbers inside a region of 
constant density diverges. The decoupling of length 
scales implies that the xy-dependence of the single 
particle wavefunctions in Eq. (12) separates in the 
thermodynamic limit. Using scaled variables, Eq. (3), 
this yields 

'Id 2 1_ 2 2gt _ mf . ' 
2dz A 2 V27rw 2 

Mf,z) = e -4^-Mf,z), (14) 

for the eigenfunctions, 0„(f, z), and eigenvalues, e^(r), 
in the confined direction, at a given radial distance, f. 
The reduced local density h™J is given by the normalized 
wavefunctions <f>(f, z): 

hf(r) = - log [1 - exp ({1(f) - e„(f)/t)] ■ (15) 

The position-dependence in Eq. (14) and Eq. (15) only 
enters parametrically through the f-dependence of the 
chemical potential, 

f 2 

M( ?, )=A-y ; (16) 

and the local-density approximation becomes exact in 
the quasi- two-dimensional thermodynamic limit. Within 
LDA, density profiles (as in Fig. 1) are directly related 
to the equation of state n(p.) of a quasi-two-dimensional 
system, which is homogeneous in the xy-planc. 

The Schrodinger equation of Eq. (14) is conveniently 
written in the basis {i/jqi V'li • • • i V'm • • • } °f t ne one- 
dimensional harmonic oscillator with u> = m = 1, as 



it diagonalizes Eq. (14) for g — 0. Using 

<M 5 ) = ^a^ip^iz), (17) 

(where we have dropped the index corresponding to f or, 
equivalently, to fi), we can write it as a matrix equation 

(A-|)a, = (18) 

with eigenvalues e u /Co z and eigenvectors a„ = 
{a 0l/ , . . . , a„„}, of the (n + 1) x (n + 1) matrix 

A Mt/ = vSpu + / —_ / dz ipfj,(z)nfd{fi,z)tpu(z) (19) 
VZiroJz J 

(where < fi, v < n) and the density n™J(jl, z) = 

— Si/ l°g[l — exp(/i — e v /t)]. The wavefunctions 
tj) v are easily programmed (see, e.g., [18] sect. 3.1), and 
the self-consistent mean-field solutions at each value of p, 
can be found via iterated matrix diagonalization. 

This full solution of the LDA mean-field equations is 
analogous to the one in Ref. [10]. The mean-field ver- 
sion used in [11], however, neglects the off-diagonal cou- 
plings in A vyL with v ^ fi. In [6], we used a simplified 
mean-field potential in order to reach explicit analytical 
expressions. These different mean-field approximations 
essentially coincide at all relevant temperatures [11], but 
ground-state occupations in z slightly differ. Further re- 
placing the coupling constant g by g tanh 1 ^ 2 (uj z /2t) has 
allowed us, in Ref. [6], to improve the agreement with 
the QMC results close to the transition. This is because 
mean-field theory overestimates the effect of the inter- 
actions in the fluctuation regime. In the following we 
always quantify beyond-mean-field corrections with re- 
spect to the full LDA solution of Eq. (19). 

III. CORRELATION DENSITY AND 
UNIVERSALITY 

Comparisons between the QMC and the mean-field 
density profiles are shown in Fig. 1 for the ENS parame- 
ters at reduced temperature t = 0.71, slightly above the 
Kosterlitz-Thouless temperature. Finite-size effects as 
well as deviations from mean-field theory are visible for 
nA T > 5. In this section, we concentrate on correlation 
corrections to mean-field theory in the thermodynamic 
limit and postpone the discussion of finite-size effects to 
Section IV. We analyze the QMC density profiles within 
the validity of the LDA, Eq. (16), and compare QMC and 
mean-field densities at the same local chemical potential 
[29] which defines the correlation density, An, 

An(/i) = n{p) - n m t{p). (20) 

As in experiments, the chemical potential is not a control 
parameter of the QMC calculation, but it can be obtained 
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from a fit of the wings of the density profile with h < i 
to the mean-field equation of state. 

Mean-field effects take into account the dominant in- 
teraction effects which, in particular, determine shape 
and energies of the ground and excited states in the 
tightly confined direction. One expects that correla- 
tion effects do not modify these high-energy modes, but 
merely affects the low-energy distribution of xy-modcs in- 
side the confining ground state, (f>o(z). This assumption 
is supported by a direct comparison of the normalized 
density distribution in z between the QMC solution and 
the mean-field approximation (inset of Fig. 1) at differ- 
ent radial distances, r. For small f, the ground state of 
the confining potential is strongly populated. For larger 
r, higher modes of the one-dimensional harmonic oscilla- 
tor are thermally occupied, and the density distribution 
broadens. However, the normalized density profile in z 
is everywhere well described by mean-field theory, and 
correlation effects hardly modify the mode structure in 
the confined direction. 



s 




position r = rfi 



FIG. 1: QMC density profile for the ENS parameters (g = 
0.13, Cb z = 0.55) at temperature t = T/tH c = 0.71 for 
different values of the particle number, compared to the LDA 
mean-field solution of Section II B at the same total number 
of particles. The inset shows the density distribution n(z) at 
different values of the radial distance f. 

In the homogeneous two-dimensional gas, corrections 
to mean-field theory at small g are described by classical- 
field theory, and correlation effects in the density pro- 
file can be expressed in terms of a universal function of 
/3(fj, — 2gn m {) [8, 12]. In a quasi- two-dimensional geom- 
etry, the corresponding relevant quantity is given by the 
local mean-field gap A m f between the local ground-state 
energy in the confining potential and the local chemical 
potential, 



A mf (f) = e (f)/t-/2(r). 



(21) 



Within mean-field theory, it fixes the local xy-density in 
the ground state of the confining potential 



log 1 



In the strictly two-dimensional limit, we have 

Amf -> (3(2gh 2 n m [/ m - //) = gn mt / ir - fi, (23) 

where n m f is the total mean-field density. Deviations are 
noticeable for large densities, as illustrated in Fig. 2 (we 
have absorbed the zero-point energy uj z /(2i) in the chem- 
ical potential). Within LDA, we expect that the correla- 
tion density An coincides to leading order in g with the 
classical-ficld-theory results of the homogeneous strictly 
two-dimensional system [12], expressed as functions of 
the mean-field gap A m f [30]. 
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FIG. 2: Mean-field equation of state and gap A m f for ENS 
parameters g — 0.13, Cb z = 0.55, t — 0.71. The gap differs 
from the approximation jl — gn/n only at high density. The 
2d mean-field curve (see Eqs (22) and (23)) illustrates the de- 
pendence of the equation of state on microscopic parameters. 

In Ref. [12], the critical density n c and chemical po- 
tential fl c at the Kosterlitz-Thouless transition were de- 
termined to 



fir 



1 ?™ r- 

log— , 
9 

9 , Cm 

-iog-e : 

7r g 



380 ±3, 



13.2 ±0.4, 



(24) 
(25) 



and the equation of state in the neighborhood of the tran- 
sition was written as 

n-n c = 27rA(X),with X = (pi- pL c )/g. (26) 

The function A(X) was tabulated. Consistent with the 
classical-field approximation, we can expand Eq. (22) to 
leading order in A m f, 



nf(A ml ) = - log(A mf ) 
so that we can express X through A m f: 



(27) 



X(A mf ) = 2=£ - 

7T g 



(22) 



_ ^£ 
9 

A m£ , 1 



9 



7T 



log i 



(28) 
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Thus, we obtain the correlation density as a function of 
A mf : 



An = 27rA[X(A mf )] + log (£„A mf /$) , 



(29) 



and a straightforward inversion of Eq. (28) allows us to 
translate the data of [12] in order to obtain the correla- 
tion density as a function of the mean- field gap. An em- 
pirical interpolation of the numerical data with < 10% 
error inside the fluctuation regime, A^ f < A m f < A 
(A£j f and A^ f are defined below), is given by 



An(A mf /ff) 



-1 + 



1 



1 



A m{ /gJ 1 + ttA^ 



,72 



(30) 



where positivity is imposed since, within classical field 
theory, the correlation density must be positive and of 
order (A mf /£T 2 for A m[ /g -» oo. 
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(0.82,0.57,0.07) 
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FIG. 3: Correlation density An vs. rescaled mean- field gap 
Amf/fl- QMC data for various interaction strengths and con- 
finements il) z are compared with the interpolation Eq. (30) of 
classical-field results [12]. 

In Fig. 3, we plot the classical-field results for the 
correlation density as a function of the rescaled mean- 
field gap. We also indicate the onset of the Kosterlitz- 
Thouless transition (u = fi c , X = in Eq. (29)) at 
A c mi /g = 0.0623, which yields Ah c = 3.164. The cor- 
relation density in the normal phase is thus finite for all 
interactions, whereas the mean- field density diverges as 
h m f = — log A^ f oc — logg for small interactions at Tkt- 
In Fig. 3, we furthermore compare the classical-field data 
for the correlation density with the results of QMC simu- 
lations of quasi-two-dimensional trapped Bosc gases with 
different coupling constants g and confinement strengths 
Cj z . The QMC data illustrates that the external trapping 
and the quasi-two-dimensional geometry preserve univer- 
sality in the experimental parameter regime. However, 
the finite coupling constant g introduces small deviations 
due to quantum corrections. 

From Fig. 3, we further see that the correlation den- 
sity is reduced to roughly 10% of its critical value for 



mean- field gaps A^ lf ~ g/^- Thus, only densities with 
^ ~ K n mf(A m f « g/ir) are significantly affected by 
correlations, and hf can be considered as the boundary 
of the fluctuation regime. In fact, perturbation theory 
fails inside this regime. For a strictly two-dimensional 
system, we have 



hf « log(7r/<?), 



(31) 



and the fluctuation regime is reached for densities h > 
hf. Outside the fluctuation regime, h < hf, mean-field 
theory is rather accurate, and can be improved pertur- 
batively, if necessary. 

To understand this criterion, which is important for 
the Kosterlitz-Thoulcss to Bosc-Einstein cross-over at 
small N (sec Section IV), we briefly analyze the per- 
turbative structure of the two-dimensional single-particle 
Green's function beyond mean-field theory [8]. Within 
classical-field theory, second-order diagrams are ultravi- 
olet convergent. Each additional higher order brings in a 
factor H 2 g/m for the interaction vertex, one integration 
over two-momenta, a factor T, and two Green's functions 
(the internal lines). Dimensional analysis of the integrals 
involved shows that each vertex insertion adds a factor 
<?/A m f. This implies that perturbation theory fails for 
g/A mf > 1. 

For lower densities, 1 < h < hf, the gas is quantum 
degenerate, yet it is accurately described by mean-field 
theory. In contrast to fully three-dimensional gases, the 
quantum-degenerate regime can be rather broad in two 
dimensions for gases with g< 1. In this regime, the den- 
sity, yet normal, is no longer given by a thermal Gaussian 
distribution. This was observed in the NIST experiment 
[2] where hf ~ 5. In Fig. 4, we illustrate this effect via 
the approximations for the tails of the distribution h(f) 
with between one and five Gaussians, as 



7T 



k„ 



= 2^ knk exp 

fc=l 



kr 2 



(32) 



where iTk is determined by the formal expansion of the 
logarithm in Eq. (15), but also appears as a cycle weight 
in the path-integral representation of the bosonic den- 
sity matrix where they can be measured (see the inset 
of Fig. 4). The successive approximations have no free 
parameters. 

Figure 5 summarizes the density profiles of a strictly 
two-dimensional Bose gas in the limit g — > where the 
classical-field calculations determine the correlation den- 
sity. At the critical temperature Tkt, the density in 
the center of the trap is critical, h(0) = h c . Correla- 
tion effects are important only in the fluctuation regime, 
f < 1.25-\/g, where h(f) > hf. However, the distribution 
of the correlation density introduces no further qualita- 
tive features to the mean- field component. The density 
profile may be integrated using the interpolation formula 
for the correlation density, Eq. (30). The critical tem- 
perature of the strictly two-dimensional Bose g (IS CIS Ob 
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FIG. 4: QMC density profile for the NIST parameters t = 
T/Tgi c = 0.75, g = 0.02, and ti 2 = 0.5, for different JV. The 
data are compared to the expansion of Eq. (32), considering 
the largest (k = 1) and the four largest terms (k = 1, . . . , 4). 
The inset shows the QMC cycle weights 7r fe for N = 100 000 
(see [6]). 



FIG. 5: Schematic density profile of a strictly two-dimensional 
trapped Bose gas at Tkt for g — ► 0. For r 2> y/g, n coincides 
with the ideal gas at TJec (inset, the classical Boltzmann dis- 
tribution e~ r2//2 is given for comparison). In the fluctuation 
regime, for r < 1.25\/g, mean-field and correlation effects be- 
come important. The density diverges as ~ log(l/<?), yet the 
correlation contribution An remains finite. 



function of the total number of particles is given by 



T, 



KT 



rp2A 

J BEC 



3 ffl„„2 9 



1 + 



Gg 



16 16tt 2 



15 



■°4 



-1/2 



This expression includes correction of order g log g com- 
pared to the mean-field estimate of Refs [6, 8]. Since 
corrections beyond classical-field theory are rather small 
(see Fig. 3), the Kosterlitz-Thouless temperature of the 
strictly two-dimensional trapped Bose gas is accurately 
described by this equation even for large coupling con- 
stants. 

For general quasi- two-dimensional gases, Fig. 5 re- 
mains qualitatively correct, but the LDA mean- field den- 
sity profile in the quasi-two-dimensional geometry must 
be used. Numerical integration of this g — > den- 
sity profile for the ENS parameters with ui z = 0.55 
and g = 0.13 leads Tkt — 0.71 T|| c , in close agree- 
ment with T KT ^ 0.70 T|| c determined in Ref. [5] di- 
rectly from QMC calculations using finite-size extrapola- 
tions. The quasi-two-dimensional transition temperature 
is smaller than the one of the strictly two-dimensional gas 
(T^ ~ 0.86 T|| c for g = 0.13). For the NIST param- 
eters (g = 0.02, Cj z = 0.5), we have T KT ^ 0.74T|| C 
from the integration of the g — > density, and QMC 
data indicate a transition slightly below this value. 



IV. FINITE-SIZE EFFECTS AND 
BOSE— EINSTEIN CROSS-OVER 

A. Central coherence 

In the normal phase, the off-diagonal elements of 
the single-particle density matrix remain short-ranged. 



so that they can be described locally. From the self- 
consistent eigenfunctions of the mean-field Schrodingcr 
equation, Eq. (12) and Eq. (13), we also obtain the off- 
diagonal reduced single-body density matrix: 



3 



exp (fl-/3ej) - 1' 



(33) 



In the local-density approximation, we can separate the 
contributions of the different transverse modes, and we 
obtain 



-Wf- - 



') 



with 



41(r;r') 



d 2 k 



A 2 e ik.(r-r') 



(2?r) 2 e / 3fi2fe2 /2m+A mf (f ; ) - 1' 



(34) 



(35) 



Here we have used that within the LDA, the density re- 
mains constant on the scale Ax, so that the mean-field 
gaps at f and r' are the same. 

At low densities, where the mean-field gap is large, 
A m f 3> 1, we can expand the Bose function in Eq. (35) 
in powers of exp (— A m f), and off-diagonal matrix ele- 
ments rapidly vanish for distances larger than the ther- 
mal wavelength At- At higher densities, in the quantum- 
degenerate regime, A m f <C 1, many Gaussians con- 
tribute, and coherence is maintained over larger dis- 
tances. In the limit A m f — > 0, we can expand the de- 
nominator in Eq. (35), exp [/?fi 2 fc 2 /2m + A mt ] — 1 « 
(3h 2 k 2 /2m + A m f, and the off-diagonal density matrix 
decays exponentially In this regime, the local mean-field 
coherence length is given by £ m f = At/V^LtA^. 
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In Fig. 6 and Fig. 7 we compare the normalized off- 
diagonal coherence function in the center of the trap 



Jdz n£>(r,z;0,0) 



c(r) 



/ dz n^(0,z;0,0) 



(36) 



from QMC calculations with LDA for the ENS and NIST 
conditions. We see that for h < hf, as in the case of the 
density profile, mean-field theory accurately describes the 
single-particle coherence. However, it is evident that at 
higher densities, n > hf, where correlation effects for the 
diagonal elements of the density matrix are important, 
mean-field theory also fails to describe the off-diagonal 
matrix elements. 

To characterize the decay of the off-diagonal density 
matrix in the fluctuation regime, n > hf, we consider 
a simple one-parameter model which neglects the mo- 
mentum dependence of the self-energies in the ground 
state of the confining potential. The single parameter of 
the model, the effective local gap A(f), is chosen such 
that it reproduces the local density of the QMC data. 
The density matrix of this "gap"-model, h^\f;r') = 

J2v ^Ai/( r > r ')4 > u(z)4>v(z'), is a straightforward general- 
ization of mean-field theory, where in Eq. (35), we replace 



A„ = 



A for v = 
A m f otherwise 



(37) 



To fix the gap A of this model, we require that the diag- 
onal elements of the density matrix reproduces the exact 
density 



Outside the fluctuation regime the gap model reduces 
to the mean-field limit. Inside the fluctuation regime, 
where a direct comparison of the coherence with mean- 
field theory is not very useful, the gap model provides 
the basis to quantify off-diagonal correlations. It cannot 
describe the build-up of quasi-long-range order at the 
Kostcrlitz-Thouless transition, but its correlation length 
£a = Ar/v^7rA(f) > £ m f bounds from below the true 
correlation length in the normal phase. In Fig. 6, we 
show that the gap model accounts for the increase of 
the coherence length inside the fluctuation regime, h > 
hf for the ENS parameters. For smaller interactions, as 
in the NIST experiment, finite-size effects qualitatively 
change the off-diagonal elements of the density matrix 
(see Fig. 7). 
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FIG. 6: Off-diagonal coherence c(r) for ENS parameters with 
t = 0.71 (main graph, h > hf) and t — 0.769 (inset, h < hf) 
compared to the mean-field prediction c m f(r) and the gap 
model of Eq. (37). In the fluctuation regime, finite-size effects 
for off-diagonal correlations are more pronounced than for the 
density (see Fig. 1). 



FIG. 7: Off-diagonal coherence c(r) for NIST parameters with 
t — 0.74 (main graph) and t = 0.769 (inset) in comparison 
with the mean-field prediction, c m f(r), and the gap model, 
c A (r), defined in Eq. (37). At t = 0.769, the total central 
density is n(0) ~ 5.1 < hf, and the system is outside the 
fluctuation regime. At t = 0.74, n(0) ~ 10.5 > hf, and the 
system is close to the Kosterlitz-Thouless transition, A m f/<? ~ 
0.08. Strong finite-size effects are evident in the fluctuation 
regime. 



B. Density profile 

Finite-size effects in the density profile are less dra- 
matic than for the coherence (see Fig. 1). Within mean- 
field theory, we have compared the density profiles of the 
finite system directly with those in the thermodynamic 
limit (LDA) , using the finite N solution obtained by the 
adapted QMC calculation described in Section II B. The 
mean-field analysis indicates that correlation effects are 
at the origin of the size-effects of the full QMC den- 
sity profiles in Fig. 1, in particular, at small system size, 
N = 1000. 
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C. Bose— Einstein cross-over 

The finite-size effects in the coherence reflect the un- 
derlying discrete mode structure of level spacing ~ hu. 
Off-diagonal properties for £a ^ K ar c cut off by the 
extension of the unperturbed ground-state wavefunc- 
tion, l r = (muj/h)^ 1 / 2 , and resemble those of a Bosc- 
condensed system with a significant ground-state occu- 
pation. Whereas in the thermodynamic limit, the in- 
teracting quasi-two-dimensional trapped Bose gas under- 
goes a Kosterlitz-Thouless phase transition, the cross- 
over to Bose-Einstein condensation sets in when A f=a 
(3hiu. If this happens outside the fluctuation regime, 
A mf < A m f < fitvjj, the Bose condensation will essen- 
tially have mean-field character. Since the temperature 
scale is given by fkv/T^ c = n/\^6N, the discrete level 
spacing is important for small system sizes, N < iVf s , 
with 



1 7T 

JV fe (A mf ) ss - 



6gH 2 (A mi /~ g y 



(39) 



For small g, close to where A m f is of order g, these 

finite-size effects trigger Bose-Einstein condensation for 
small N. In particular, for systems with TV < Nf s (A^ n{ ) « 
ir 4 /(6g 2 ), a cross-over to a mean-ficld-likc Bose conden- 
sation occurs[31], whereas for N > iV fs (A^ f ) f=a 400.g -2 
Kosterlitz-Thouless-like behavior sets in (see inset of 
Fig. 10). We notice that the finite-size scale Nf s oc l/g 2 
diverges very rapidly with vanishing interactions, which 
could make the cross-over experimentally observable. 

For a finite system with N < iVf^A^j), the condensate 
wavefunction does not develop immediately a Thomas- 
Fermi shape, but remains close to the Gaussian ground- 
state wavefunction of the ideal gas with typical extension 
l r = (rriLO /h)^ 1 / 2 . Thus, for small condensate fraction 
no, deviations of the moment of inertia / of the trapped 
gas from its classical value, I c \ = J d 2 r r 2 n(r) ~ Nl 2 -,, 
are negligible, of order (I c i — I)/I c i ~ n l 2 /l^ ~ N^ 1 / 2 . 
Only for larger condensates with gNo 3> 2ir, the self- 
interaction energy dominates the kinetic energy, and the 
condensate wavefunction approaches the Thomas-Fermi 
distribution of radius ~ It, resulting in a non-classical 
value of the moment of inertia. In this low-temperature 
regime, the system can be described by a condensate with 
a temperature-dependent fluctuating phase [13]. There- 
fore, for small systems, a non-classical moment of iner- 
tia only occurs at lower temperatures than condensation, 
roughly, at a condensate fraction no > g. 

To illustrate the cross-over between the Bose-Einstein 
regime at small N and the Kosterlitz-Thouless regime at 
large N, we have calculated the condensate fraction and 
condensate wavefunction for the ENS parameter in Fig. 8. 
To determine both quantities in inhomogeneous systems, 



l 3d 



must be explicitly diagonalized, as the eigen- 



functions of the single-particle density matrix are not 
fixed by symmetry alone. In quasi-two-dimensional sys- 
tems, the full resolution of the off-diagonal density ma- 



trix in the tightly confined z-direction is difficult. It is 
more appropriate to consider the in-plane density matrix, 
rv- 1 ' (r, r'), where the confined direction is integrated over. 
?i( 1 )(r, r') = J dz J dz' n^J(j, z; r', z'). Because of rota- 
tional symmetry, n^^(r, r') is block-diagonal in angular- 
momentum Fourier components 



,.(!) 



(r,r') 



n=0 1= 



£ N nl ^ nl {r') Vnl {r)e 1 ^ 



(40) 



where a(r,r') denotes the angle between r and r', and 
N n i is the occupation number of the normalized eigen- 
mode tp„i. The (in-plane) condensate fraction, hq = 
Nqq/N, corresponds to the largest eigenvalue with I = 
and condensate wavefunction <fio(r) = </?oo( r )- Projection 
on the Fourier components is convenient for determining 
the condensate fraction and wavefunction within QMC. 

For the ENS parameters in Fig. 8, the central density 
is already inside the fluctuation regime, and the conden- 
sate wavefunction differs from the Gaussian ground state 
of an ideal gas. However, for small systems, N < 10 3 , 
it still has a Gaussian shape indicating that the conden- 
sate kinetic energy dominates the potential energy. The 
condensate fraction vanishes as uq ~ N~ x / 2 . 

The QMC calculations of [5] demonstrated that the 
condensate fraction of the quasi-two-dimensional Bose 
gas vanishes in the normal and in the superfluid phase 
in the thermodynamic limit, N — > oo. However, in the 
low-tempcraturc superfluid phase, the condensate frac- 
tion approaches zero very slowly with increasing system 
size, so that an extensive condensate remains for practi- 
cally all mesoscopic systems. 
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FIG. 8: Condensate fraction no (inset) and wavefunction 
ipo(r) (main graph), computed by QMC for ENS parameters 
with t — 0.71 for various system sizes. The condensate wave- 
function is modified with respect to the ground-state wave- 
function tpo(r) of the unperturbed harmonic oscillator, but 
the Gaussian shape is preserved. The two largest systems are 
above the scale iVf s (see Eq. (39)) 
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V. TWO-PARTICLE CORRELATIONS 
A. Pair-correlation function 

Density-density correlations can be analyzed by con- 
sidering the three-dimensional pair-correlation function, 
rv- 2 '{r ; r '). This quantity factorizes within mean-field 
theory into terms described by the one-particle density 

matrix, n| i j(ri;r 2 ) (see Section IV A): 



n 



(2) 



'{(ri;r 2 ) = n mi (r 1 )n mi (r 2 ) + n^ { (fi;f 2 ) ■ (41) 



For vanishing distances 
r , r ) 



diets n^( 



2n 2 mf (r). 



-> r 2, mean- field theory pre- 
For Bose-condensed atoms, 
this bunching effect is absent. 

In two dimensions, deviations from 2n 2 of the pair- 
correlation function at contact signal beyond-mean-field 
fluctuations [12, 19, 21]. In Ref. [12], the universal char- 
acter of the contact value was used to define the quasi- 



condensate density, n qc (r) = [2[n(r)] : 



>(r,r)] 



1/2 

V 7 ' /J 

This quantity has been studied in Ref. [22] for a quasi- 
two-dimensional trapped gas within classical-field theory. 

The pair-correlation function of the quasi-two- 
dimensional gas is obtained by integrating both coordi- 
nates over the confined direction 



,(2) 



( r i; r 2) 



dzi I dz 2 n (2) (ri,zi;r 2 ,2 2 ), (42) 



and mean-field expressions for this quantity follow from 
Eq. (41) together with Eqs (34) and (35). Inside the 
fluctuation regime, the gap model (using Eq. (37) in the 
mean-field expressions) again leads to an improved pair- 
correlation function . 

Figure 9 illustrates that outside the fluctuation regime 
mean-field theory describes the pair-correlation function 
well. In contrast to a strictly two-dimensional gas, the 
contact value of the pair correlation function is below 2. 
Even in the mean-field regime, the occupation of more 
than one mode in the confined direction causes a notice- 
able reduction of the pair-correlation function at contact. 
The above definition of the quasi-condensate must there- 
fore be modified in this geometry to maintain its univer- 
sal character. 

At short distances r ~ ro, the pair-correlation function 
depends on the specific form of the interaction. This 
cannot be reproduced by the single-particle mean-field 
approximation. However, two-particle scattering proper- 
ties dominate for small enough distances as, for example, 
the wavefunction of hard spheres must vanish for overlap- 
ping particles. This feature can be included in mean-field 
theory by multiplying its pair-correlation functions by a 
short-range term X2d(i r ), which accounts for two-particle 
scattering [23]. In two dimensions, \2d shows a charac- 
teristic logarithmic behavior for short distances: 
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FIG. 9: Central pair correlations, g(r) = n' 2 ' (r, 0)/[n(r)n(0)], 
of the quasi- two-dimensional trapped Bose gas at temperature 
T = Tjlc (main figure) and T = 0.75 T|| c (inset) for N = 
100, 000 atoms (ENS parameters) , together with the predic- 
tion of the mean-field gap model, gA{r) = (r, 0) /n(r)n(0), 
and the short-range improved mean-field model, X2d(f)gA(r). 



Factorizing out the short-range behavior from the pair- 
correlation function, the correlation part of the renor- 



malized pair-correlation function, h A 



/X2d, should 



be dominated by contributions from classical-field the- 
ory, its contact value is universal, and it might be used 
to define a quasi-condensate density in quasi-two dimen- 



sions via n 2 

4 U 



lim r . 



"a } M) 



n (2) (r,0)/xM(r) 
(see Fig. 10). Similar to the correlation density, the 
quasi-condensate density is universal. At Tkt, the clas- 
sical field result is n qc ~ 7.2, whereas it is around 2.7 at 
the onset of the fluctuation regime, so that, in the normal 
phase, n qc /n vanishes as | logg| _1 for g — > 0. 



B. Local-density correlator 

Due to the three-dimensional nature of the underly- 
ing interaction potential, observables which couple di- 
rectly to local three-dimensional density fluctuations in- 
volve the following density correlator 



K 



(2) 



lim 

5^0 



2nl z J d2 



,,(2) 



(r, z; r + 5, z) 



X3d(S) 



(43) 



where X3d(r <C At) — (1 — a s /r) 2 describes the uni- 
versal short-distance behavior of the three-dimensional 
two-body wavefunction in terms of the s-wave scattering 
length a s . 

Arguments similar to those in Section VA show 
that the local-density correlator in general differs from 
the contact value of the quasi-two-dimensional pair- 
correlation function, Further, the integration over the 
square of the ground state density in z leads to a 
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FIG. 10: Quasi-condensate density n qc obtained by QMC 
from the renormalized pair-correlation function at T = 
0.71 Tii c , and T = 0.75 T|| c (see Fig. 9) for ENS parame- 
ters, plotted as a function of A m f, and compared to classical- 
field simulations [12]. The inset shows the boundary of the 
region with strong finite-size effects (see Eq. (39)). 

phase-space-density dependence which destroys the sim- 
ple mean-field property oc 2n of strictly two- 
dimensional Bose gases. 

VI. CONCLUSIONS 

In this paper, we have studied the quasi-two- 
dimensional trapped Bose gas in the normal phase above 
the Kosterlitz-Thouless temperature for small interac- 
tions g < 1. We have discussed the three qualitatively 
distinct regimes of this gas: For phase-space densities 
h < 1, it is classical. At higher density, 1 <n<fif, the 
gas is in the quantum mean-field regime, and its coher- 
ence can be maintained over distances much larger than 
At- Finally, mean- field theory fails in the fluctuation 
regime hf < n < h c and beyond-mean-field corrections 
must be taken into account. 

In the fluctuation regime, for small interactions, the 
deviations of the density profile with respect to the mean- 
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field profile are universal. Mean-field theory thus ac- 
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